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NOMENCLATURE 


Flexible and rigid areas of cavity 
Plate length 
Plate width 
Cavity sound speed 

Eh3/ 12fl-\' 2 ) 


Modulus of elasticity 
th 

in cavity natural mode 
Plate admittance function 
Plate thickness 

Plate transfer function; see Eq. (28) 

= jn w 2 a 4 
n 

D 


see equation (62), et. seq. 

Characteristic lengths of random pressure field 
Plate generalized mass; also cavity generalized mass 
Plate mass/area 
normal 


= pa 4 ; see equation (57) 
hD 


Pressure on plate 
Generalized force on plate 
Generalized plate coordinate 
Correlation function 
Time 

Pressure field convection velocity 
Cavity volume 
see equation (57) 



Ill 


w 

x, y, z 
a, 6, Y 


V 2 

5 * 


m 


m 


m 


Plate deflection 
Cartesian coordinates 

Empirical constants for random pressure field 
Laplacian 

Boundary layer thickness 

Power spectral density 

Plate density 

stress 

Dummy time 

Modal damping 

Modal frequency 
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1 . INTRODUCTION 

In the present report we shall treat the response of a structure 

to a convect ing-decaying random pressure field. The treatment 

follows along conventional lines after Powell and others. That is, 

the pressure field is modelled as a random, stationary process 

whose correlation function (and/or power spectra) is determined from 

experimental measurements. Using this empirical description of the 

random pressure, the response of the structure is determined using 

[2 31 

standard methods from the theory of linear random processes. 5 
The major purpose of the report is to provide a complete and 
detailed account of this theory which is widely used in practice 
(in one or another of its many variants). A second purpose is to 
consider systematic simplifications to the complete theory. The 
theory presented here is most useful for obtaining analytical 
results such as scaling laws or even, with enough simplifying 
assumptions, explicit analytical formulae for structural response. 

Some of these latter results are thought to be new; however, so 
many authors have considered various simplified versions of the 
general theory, the authors hesitate to claim novelty for any 
specific result. Hence, the title of the present report. 

It should be emphasized that, if for a particular application 
the simplifying assumptions which lead to analytical results must be 
abandoned, numerical simulation of structural response time histories may 

[4 5] 

be the method of choice. 1 1 J Once one is committed to any sub- 
stantial amount of numerical work (e.g. computer work) then the 
standard power spectral approach loses much of its attraction. 
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2. INPUT-OUTPUT RELATIONS FOR STRUCTURE 

In the present section we derive the input-output relations for 
a flat plate. It will be clear, however, that such relations may be 
derived in a similar manner for any linear system. 

The equation of motion for the small (linear) deformation of a 
uniform isotropic flat plate is 

DV 4 w + m a 2 w = p (1) 

Tt? 

where w is the plate deflection, p the pressure loading and the 
other terms are defined in the Nomenclature. Associated with (1) 
are the natural modes and frequencies of the plate which satisfy 

DV 4 ^ - w n 2 m i^ n = 0 (2) 

where w is the frequency and (x,y) the shape of the n th 
n n 

natural mode. In standard texts it is shown that the i^ n 
satisfy an orthogonality condition 

// dxd y = 0 ^ 

for m ^ n 

If we expand the plate deflection in terms of the natural modes 

w = Z q n (t) ^ n (x,y) (4) 

n 

then substituting (4) into (1), multiplying by i|> m and 
integrating over the plate area we obtain 



3 


M m [q m 


m 


“J - <3m 


m = 1,2, • •• 


(5) 


where we have used (2) and (3) to simplify the result, 
are defined as 


M and 
m 


M = / / m il 2 dxdy 
m ; J m 

<*m s ffP+ n dxd >' 


( 6 ) 



For structures other than a plate the final result would be un- 
changed, (5) and (6); however, the natural modes and frequencies 
would be obtained by the appropriate equation for the particular 
structure rather than (1) or (2). Hence, the subsequent 
development, which depends upon (5) only, is quite general. 

Before proceeding further we must consider the question of 
[structural) damping. We shall defer a discussion of acoustic or 
fluid damping to a subsequent section.- Restricting ourselves to 
structural damping only we shall include its effect in a gross 
way by modifying (5) to read 


M [q + 2c;w q + cu 2 q ] = (7) 

where is a (non-dimensional) damping coefficient usually 

determined experimentally. This is by no means the most general 
form of damping possible. However, given the uncertainty in our 
knowledge of damping from a fundamental theoretical viewpoint 
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(see [6]) it is generally sufficient to express our meager knowledge. 
If damping is inherent in the material properties (stress-strain law) 
of the structure, the theory of viscoelasticity may be useful for 
estimating the amount and nature of the damping. However, often 
the damping is dominated by friction at joints, etc., which is 
virtually impossible to estimate in any rational way. 

Now let us turn to the principal aim of this section, the 
stochastic relations between loading and response. We shall obtain 
such results in terms of correlation functions and power spectra. 

The correlation function of the plate deflection w is defined 
as 

T 

R (i;x,y) = lim 1_ / w (x,y, t)w (x,y, t+x)dt (8) 

2T -T 

Using (4) we obtain 


R w 0 ; x, y) 


E E i^ m (x,y)^ n (x,y) R (t) 
m n 4 m H n 


(9) 


where 



(t) = lim 
T->-“ 


T 


i_ / 

2T -T 


q m (t)q n (t+T)dT 


(10) 


is defined to be the cross-correlation of the generalized coordinates, 

v 


Defining power spectra 
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V w ;x,y) = I 

TT 


$ (u>) = 1 

q q — 

n m n n tt 


/ R w (t;x,y) e ' dx 


/ R (x) e dr 

; q q 
-co n m n n 


( 11 ) 


( 12 ) 


we may obtain from (9) via a Fourier Transform 


* w (^;x,y) = l l 
m n 


4 m H n 


(13) 


(9) and (13) relate the physical deflection w to the generalized 
coordinates or displacements q^. 

Consider next similar relations between physical load p 
and generalized force Define the cross-correlation 

T 



lim 1__ / 0 (t)0 (t+T) dt 

2T -T 


(14) 


Using the definition of generalized force (6) 


Q^t) = // p(x,y,t) ^ m ( x >y) dxd >’ 

Q^(t+x) = // p(x*,y*,t+i) ip n (x*,y*) dx*dy* 


and substituting into (14) we obtain 
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R Q m Q (t) s //// \(x,y)^ n Cx*,y*) 


Rp(x;x,y,x*,y*) dxdydx*dy* 


where we define the pressure correlation 


(15) 


R (t ; x, y, x* ,y*) = lim 1_ / p(x,y, t)p(x*,y*, t+i)dt 
P 2T -T 


(16) 


Note that a rather extensive knowledge of the spatial distribution 
of the pressure is implied by (16). 

Again defining power spectra 


Vo M - i ! Vo CT) elMT dT 

in ti 7T -« in ti 


(17) 


* (“jx, y,x*,y*) = 1_ J R (T;x,y,x*,y*) e U)T dr 


TT -<» 


we may obtain from (15) 


\Qn (a ° = //// ^ m (x,y)^ n (x*,y*) 


(18) 


(19) 


$p (u>;x,y,x*,y*) dxdy dx*,dy* 


Finally, we must relate the generalized coordinates to the 
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generalized forces. From (7) we may formally solve (see [2]), for 
example) 


oO 

= ! H n W dt l 

.CO 


( 20 ) 


where 


CO 

H (t) = 1 / I (oi) e lwt d w 

2 it - OT 


( 21 ) 


and the "transfer function", 


I (ui) = M fa) ^ + 2 \ iw uj-uj^ ] 
n n L n n n 


Also 


I M « / H n (t) e' lut dt 

.00 


which is the other half of the transform pair, cf (21) 
From (20) and (10) 


R M - lim 1 /// ^Ct-tp H n (fT-t 2 )Q 11| Ct 1 }Q n (t 2 ) 

T-* 00 2T -T 


dt 1 dt 2 dt 


Performing a change of integration variables and noting (14), 
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R Vn (T) = -» HmCO HnCn3 R Q.Q n tT ‘ n+ ° d5dn 


( 22 ) 


Taking a Fourier Transform of (22) and using the definitions of 
power spectra (12) and (17), we have 


* l m M V -* 0 


(23) 


Summarizing, the relations for correlation functions are (9), 
(15), and (22) and for power spectra (13), (19) and (23). 

For example, substituting (19) into (23) and the result into 
(13) we have 


4> w (co;x,y) = 2 z if> m (x,y) ^(x.y) I m (o0 I R (-w) 


m n 


//// Cx,y) (x* ,y*) 


(24) 


(w;x,y,x*,y*) dxdydx*dy* 


This is the desired final result relating the physical loading 
to the physical response in stochastic terms. 

Sharp Resonance or Low Damping Approximation 

Often (2 A) is approximated further. Two approximations are 
particularly popular and useful. The first is the "neglect of 
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off-diagonal coupling". This means omitting all terms in the double 
sum except those for which m = n. The second is the "white noise" 
approximation which assumes that 4^ is essentially constant relative 
to the rapidly varying transfer functions I (go) . Making both of 
these approximations in (24) we may obtain the mean square response 

oo 

w 2 (x,y) e R (i=0;x,y) = / 4 fto;x,y) d to 

w o • 

//// v x >y) V x *' y * } 

•4p (w m ;x,y J x*,y*) dxdydx*dy* (25) 

Of course, only one or the other of these approximations may be 
made, rather than both. However, both stem from the same basic 
physical idea: The damping is small and hence, 1^ has a 

sharp maximum near m = tu . That is 


^ n Z 4 2 (x, y) 

m 

4 m 


M 2 ui 3 4 
m m ^m 


I (to ) I (-to ) << 
nr n^ m 


I (to ) 
m m 


2 


I (to ) I (-to ) << I I (to ) ! 2 
m v n n v n J 1 n v n' 1 

and the "neglect of off-diagonal coupling" follows. 
Also 

/ ^ I I 2 dto £ $ p U> m ) / | I m (o0 I 2 du 


and (25) follows by simple integration. 
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Note that if we take the spatial mean square of (24) then 
using orthogonality (for a uniform mass distribution) one may show 
that the off-diagonal terms do not contribute (see Powell^). 

Finally note that if we desire stress rather than deflection, 
then it may be shown that analogous to (25) one obtains 



s ° m 2 ( x >y) //// ^ m ( x *»y*) 


• (u> m ; x,y,x*,y*) dxdydx*dy* (26) 


where a is stress due to w = . 

m m 
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3. SIMPLIFIED RESPONSE CALCULATIONS 

Miles 1 Approximation for_ Spatially Well-Correlated Noise Sources 
(Acoustic source characteristic length is large compared to 
structural dimensions) 

One sometimes makes the further simplifying assumption that 
the acoustic pressure loading is perfectly correlated over the 
plate, i.e. 

$ p (w; x,y,x*,y*) = <$ p (to) (27) 

and the power spectra is independent of spatial coordinates. 

This is reasonable provided the characteristic length associated 
with the acoustic source is large compared to the plate length 
and width. For jet engine noise the characteristic length is 
on the order of the engine diameter and this assumption is par- 
ticularly useful. For boundary layer noise, where the charac- 
teristic length is much smaller (on the order of the boundary 
layer thickness), it is less so. It is normally a conservative 
assumption in that the plate response is overestimated. 

Using (27), (25) and (26) we see that 

1 = //// $p <j> m (x»y) dxd y dx*dy* = 

w H uj v* xd ^ 2 < 28 > 

Hence (25) and (26) become 
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w 2 (x,y) = tt_ £ t 2 (x,y) * (w ) F/J dxdyl 2 (29) 

4 m PL J 


M 2 co 3 C 
m m m 


^ (x,y) = it i a m 2 ( x »y) VV [F/ dxd y ] 2 

4 m PL J 


(30) 


M 2 w 3 ? 
m m m 


These may be interpreted in a particularly helpful way by 
recognizing that for a uni form , unit , static pressure load 
(see (5) and (6)) 



M u> ' 
m m 




M u) 
m m 


2 


(31) 


and thus the physical deflection due to this uniform, unit, static 
load is 


w S = E q ip = £ i p [fty dxdy 
M m y m r nr J m } 

m m rr ■ — 7 — 

M a] * 
m m 


and the corresponding stress is 


cr S = £ q a = E a //V dxdy 
n m m nr * m 

m m rj -? — 

M u z 
m m 


(32) 


(33) 


If the response is dominated by a single mode (both acoustically 
and due to the static load) , then the summations in (29 - 33) 
may be ignored and (29) and (30) written as 
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w 2 = jVjl ( V (wS)2 


m 


(34) 


= ir u 

4-^ 


4 rio ) 
J2. m 


(O 


m 


(35) 


r 7i 

The above forms are due to Miles 1 

Nearly Uncorrelated Noise Sources 

(Acoustic source characteristic length is small compared to 
structural dimensions) 

For some noise sources, e.g. turbulent boundary layers, the 
assumption of perfect correlation is not satisfactory and a 
different limiting approximation may be more useful. A typical 
(empirical) form of the pressure power spectra is^ 

- | x-x* [ - | y-y* | 

4 (w; x, y,x*, y*) = 4 n (w) e L e L cos w | x-x* | (36) 

P u A / U 

c 


We wish to consider the case where L << a, L << b and 

x y 


L ,L - (boundary layer) noise characteristic lengths in 
x y 

x, y directions 

a, b - structural lengths in x,y directions 

In Miles’ approximation, it was implicitly assumed that L >> a, 

X 


L >> b. 

y 
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If we substitute (35) into (25) or (28), the result is 


1 - //// 4 m ( x .y) ^ m (x*,y*) dxdy dx*dy* (36) 

% * 0 M 4 L x L y // ip m 2 (x,y) dxdy 

To obtain the above we note that if L << a, L << b then 

x y 

only for x* ■v x and y* ^ y will the integrand of right hand 
side of (36) contribute. Thus 

I * W//* m 2 (x,y)j//e-|^ e-|^*-l cos “I255I 

Performing the integral over x*» y* (for definiteness consider a 
rectangular plate of dimensions, a,b) 


dx*dy* I dxdy 


I v $ ( U )L L // ijj 2 j[2-e" x/ L cos kx/ T - e“ (a_x)/ L 
'v 0 ^ J x y J J Y m ) 1 x L 3 


cos k (a-x)/ L ] 
x 


[2 - e" y/ L y - e- (b - y)/ L y ] dxdy 


where 



_ x / 

Neglecting e L , etc. compared to 2, one finally obtains 
(36). Note this last step is a particularly good approximation 
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\ 

if ^ is small near x = 0, a and y = 0, b. Conversely, if 
had its maxima near the plate edges, the above might require 
modification (this is not the usual case). 

If one still assumes the response is dominated by a single mode 
(a more questionable assumption for L << a than L >> a ) 
we see, comparing the above result for I to Miles' approximation, 
that the latter may be modified to give the structural response 
by multiplying by a factor of (square root for rms response) 

4 L x L y // ^ 2 (x,y ) dxdy (37) 

[//i^dxdy] 2 

For typical this ratio will be on the order of 

4(L /a) (L /b) and hence, as expected, leads to a reduction 
x y 

in the estimated structural response from that given by the 
Miles' approximation. 


1 Miles 
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Numerical Example 

As an example, we have considered a clamped isotropic plate 
(with and without tiles) representative of those recently tested 
in a wind tunnel for Space Shuttle^. Although this was nomi- 
nally a panel flutter test, some tile failures occured for reasons 
other than flutter. Hence, it is of interest to assess the possible 
impact of vibroacoust ic response on these panels. 

Metallic Panel Only: 

The aluminum metallic panel has the following properties 


H = 10 7 psi 

a = 

24" 

P * .1#/. 3 

b = 

18” 

m m 

h = 

.05 


The first natural mode is approximated by 


= (~1 - cos 2ttx”| fl - cos 2 iry 

L a J L b 

Using a Rayleigh-Ritz approximation the calculated natural 
frequency is 37 cps as compared to the measured value of 42 cps L J . 
Also it is calculated that 

w s = .62 inches/psi 

and c s = 2360 psi/^^ 
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From Figure 12 of Reference [9], = .00305 psi 2 / c p s at 37 cps 

and it is assumed that = .01. 

Thus, using Miles' approximation. 


\! a 2 = VT [ w i ^ 0,5 


= 7400 psi 


However, from Vaicaitis ^ , et. al.. 


L = 1.22 6* , L = .26 6* 

x y 


where 6* is the boundary layer thickness. Taking 6* = 1 inch 
and using the previously assumed , we compute from (37) 


= 9 L L 
_* JL 


= .0066 


*Miles 


which reduces the rms stress to 


\/ a 2 = 7400 x [.0066] 1//2 = 600 psi 


Addition of Tiles: 


The measured value of cai is 85 cps with 1" tiles resting on 
a soft felt strain isolator which in turn is bonded to the metallic 
panel^. From Muhlstein^ , et. al. , = .00148 at 85 cps. 
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It is assumed that = .02 and 

o S .^, . is estimated from 

with tiles 


s 

with tiles 


“ Cmass x Ul 2 ) tiles ° n0 tilM 


(mass x wi 2 ).. 


tiles 


(The above could be computed more precisely from the analysis of 
Reference [10] if it were necessary.) 

Now 


mass , - .005 #/ . 2 

no tiles m 


mass^ . 1 
tiles 


» (.005 + .0052) #/. 


2 _ 


m 


metallic plate and tiles 


= 295 psi/ 


psi 


The stress in the metallic plate with tiles added is then 



655 psi (Miles* approximation) 

53 psi corrected for L , L effect. 

x y 


The bending stress in the tile itself is computed to be 65.5 psi 
and 5.3 psi respectively. Here, we make the conservative 
assumption that the isolator is ineffective and the tile bends 
with the metallic plate, i.e. 
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Tile 

a 

s 

■ t^Tile 

= 50,000 

x 1 

Aluminum 

a 

s 

Aluminum 

10 7 x 

.05 


The through the thickness stress (assuming tile moves up 
and down on isolator with no motion of metallic plate) is 
computed using the following values 

a = 1 psi/ psi 

Cl = - 1 


co i = 680 cps 


<i> = .0002 psi 2 / 

p cps 


The result is 


o 2 = 1.2 psi Miles’ Approximation 


Correcting for L , L but now using dimensions of tile, a = b = 6", 

* y . 

we have 


.079 


and 


"'‘Miles 


V* 


1.2 x t .079] 1/2 = .34 psi 
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Tile 

°s 

Aluminum 

o 

s 


(El) Tile 

= 50,000 x 1 

^ Et ^Aluminum 

10 7 x .05 


The through the thickness stress (assuming tile moves up 
and down on isolator with no motion of metallic plate) is 
computed using the following values 


= 1 psi/ 


psi 


Cl = 


a) i = 680 cps 

*p " - 0002 P si2/ ops 
The result is 



1.2 psi 


Miles’ Approximation 


Correcting for L , L but now using dimensions of tile, a = b = 6", 
x y 

we have 


079 


Miles 


1/2 _ 


and 


1.2 x t • 079] 


.34 psi. 
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Trial Modes 

Sometimes it is convenient to use assumed trial modes rather 
than the natural modes. To determine the latter may itself involve 
a rather elaborate calculation. 

Define 

ib - natural mode o - stress due to A 

m mm 

d>. - trial mode a. - stress due to 

J J 3 

where, assuming for simplicity that $ (u> m ) is independent 

of position, 




s [//% dxd /3 2 


m_ 

M 2 tu 3 s 
m m m 


(38) 


M = f/m <p 2 dxdy 

m ' ' m 

- m t ^ 1 natural frequency 

r - m^ damping ratio 
m 

Clarkson and others advocate retaining only m = 1. 
It will be shown (see equations (41) - (48) that 


/ / ijydxdy = bf a b 

M = mab Z (b 1 ?) 2 
m . J 


( 39 ) 
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(38) and (39) may then be used to calculate a 2 , viz 


a 2 = £ q 7 a. 

H m m 
m 


( 40 ) 


Here we neglect q^q^. OPTION I. 


Relationships Between Natural Modes and Trial Modes 

We have 


ib = I b. m <J>. 

m j J J 


(41) 


where bT are the eigenvectors of a natural mode calculation 
using the trial modes, <t>y 
We also know 


w = £ 
m 


% * 


m 


; a j ** 


(42) 


Using orthogonality property * of we determine from 

(42) that 


% = * a j U *m dxdy 

// ^ m 2 dxdy 


C43) 


Ising (41) in (43) and orthogonality property * of <|k 




£ a. b. 

j 3 3 

£ (b , m ) 2 

i 3 


m 


C44) 


*Note we are assuming mass distribution is uniform. 
**a. are generalized coordinates associated with 4> 
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Thus 

q q = Z Z a. a, b.^b, 11 (45) 

m n j k J k 

[E (b V][I OvV] 

J k K 

In practice we probably only need consider m = n. 

We also need for some purposes to consider 

M 5 // m ^ m 2dxd y and // V* xdy 

Using (41) 

// m 2 dxdy - m ZZ b, m b. m // <j>, <t>.dxdy 
m kj K - 1 K J 

or M = m I (b. 111 ) 2 a b (46) 

m j J 

where a,b are length, width of plate 
and // i^dxdy = Z b^. m // Q^dxdy 

= b” a b (47) 

Also from (40) and (41) 

a = Zb > m a. (48) 

m j 3 J 

cf. (39) and (46) - (48). 

* again assuming mass is uniformly distributed. 
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Summary: • 

Note this method is not limited to uniform pressure loading, 
even though for convenience tve have made this assumption up to 
this point. 

OPTION II 

Calculate 

T 

a. a, = lim 1_ / a. (t)a. (t)dt (49) 

J K x-wo 2T -T •' K 

The a. would have to be determined from a suitable modal 
J 

dynamic analysis. 

Calculate 


• Z . I °j°k 

3 k J J 


o j - stress at x, y due to w = 4 . (x,y) 


OPTION III 


Calculate a. a, 
3 k 


Calculate q 

ti 


m , n 

£ l a. a. b. b, 

i k ^ k 3 k 

tI(b. m ) 2 ][E(b n ) 2 ] 
4 J V K 


Calculate 


? = H q q a a 
n m n n m n 

m n 


m 


x, y due to w = ip (x,y) 


(50) 


stress at 
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Since are natural inodes, it should be possible to approximate 

qq 'v- q ^ 6 

mi^n = n m mn 

Note b^ m are eigenvectors from a trial mode eigenvalue 
calculation. 
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4. SOUND PRESSURE LEVEL IN AN ACOUSTIC CAVITY 


In many instances the transmission of random pressure 
fluctuations from an external flow through a flexible wall into 
an interior cavity is of concern. Given the wall motion, the 
cavity pressure (sound) field may be determined. However, there 
also may be a "feedback" in that the cavity pressure modifies 
the plate motion. There also may be "feedback" due to pressure 
changes in the external flow as a result of flexible wall motion. 
Both give rise to acoustic or aerodynamic damping which may 
substantially modify the wall motion. Although this modification 
may be significant^ in some instances, in this elementary dis- 
cussion we ignore it. 

The following analysis follows closely the derivation of 
Ventres ^ which is notable for its generality and simplicity. 

The analysis is valid for arbitrary cavity geometry and time 
dependent motion. The end result is an expansion for the acoustic 
pressure within the cavity due to motion of the walls, in terms 
of the normal modes of the cavity with its walls assumed to be 



Let the cavity occupy a volume V, and be surrounded by a wall 


surface A, of which the portion Ap is flexible, while the 


remainder 


*R 


is rigid. 


If the fluid within the cavity is at 


rest prior to motion of the wall, the fluid pressure p satisfies 




26 - 


the familiar wave 

equation, and associated boundary condition: 

v2 p - l / , i 

; c 2 8 

£ “ 0 

(51) 

3p= -p a 2 w 
an at" 2 

on A p 

(52) 

= 0 

on A r 


In these equations p and c are 

the fluid density and acoustic 

velocity within the cavity, and w 

is the displacement of the 

flexible portion 

of the wall in the normal direction n (positive 

outward) . 



Equation (51) has normal mode 

solutions F e lwt , m = 0,1,2 . 
m 

with the following properties 


V 2 F = ~ ( 

2 

U ) \ r 

— 1 m 
c / 

(53) 

aF =0 

m 

on A 

(54) 

an 



1 / F F dV 

v V m n 

= 0 m ^ n 

(55) 


- M men 

n 



men 
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Note that equation (53) has the solution co Q = 0, F Q = 1. All 
other frequencies , m = 1,2, . . . are positive, however. 

The wave equation (51) can be transformed into a set of 
ordinary differential equations in time by using Green's 
Theorem in the form: 


/ (PV 2 F - F V 2 P) dV = 
m m 



A 



F 

m 



3n 


(56) 


By defining 


p„ = / p p m dv 


m 


pc 2 V 


W = 1 J w F dA 
m V A p 


(57) 


and making use of the fact that p and F m satisfy equations (51) 
and (53), and boundary conditions (52) and (54), the following 
ordinary differential equations are obtained from (56) : 


A dot ( ) 

quantities 

expansions 


u> 2 p = - W (58) 

mm m 

denotes differentiation with respect to time. The 

P and W are the coefficients in normal mode 
m m 

for the pressure and plate deflection, e.g. 




1 P 


m 


m 


M 


m 


F 

m 


P 


(59) 
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Since the normal modes F m satisfy the homogeneous boundary con- 
dition (54) on the entire wall surface A, the normal derivative 
of expression (59) does not converge uniformly on the flexible 
portion, A p , of the wall surface. Expression (59) is suitable, 
however, for calculating the pressure itself throughout the 
cavity and everywhere on the wall surface, including the 
flexible portion. 

Equations (58) may be easily solved for arbitrary wall 
motions. The solution appropriate for the usual initial con- 
ditions, p = 3 jd = 0, at t = 0 is 

at 


p 0 - - tyt) 

- - l V T > cos “ ra (t ' T) dT 


(60) 


In keeping with the intention that the cavity be in equilibrium 
at t = 0 and for all previous time, it has been specified that 
w = 0 and 3w_ = 0 at t = 0 in (60) . 

at 

The wall deflection w is often expressed as a series of 
the form 

w = 1 Sn C61) 

m 

in which the modal functions are defined over the region 

Ap, their properties being determined by structural considerations. 
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In this situation the quantities are, from (57) 


W = A c L q 
m F mr n r 


( 62 ) 


where 


L =1 / F MA 

mr A“ A m r 
F F 


Equations (58) become 


P + w 2 P 
m mm 



L 

mr 


(63) 


In the present context the relationships between structural 
motion and cavity pressure power spectra are of interest. 

Taking a Fourier Transform of (63) , we have 


( '“ 2 % 2) P m = ^ L mr“ 2 ? r (64) 

V 

Defining power spectra in the usual way, we obtain 


4> 

= |h 

12 * 


C65). 

P P 
m n 

1 mnrs 

1 i 

V*s 

where 

H = 

mnrs 

H mr 

H 

ns 

C66) 

and 

H = 

mr 

a f 

V 

L io 2 

mr 

(-u) 2 +u m 2 ) 

(67) 


Alternatively using (60) and (62) one can relate the correlation 
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functions of P m and q r * Subsequently, one can relate q r 
to w via (61) and to P via (59). Eventually, we may obtain 
a relationship between the power spectra of P and w, the desired 
result. Since much of this analysis parallels that of earlier 
sections on the relationship of external pressure to structural 


motion we omit details here. 
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5. EMPIRICAL EXPRESSION FOR RANDOM PRESSURE SPECTRA 

From (25) it is seen that the pressure power spectra is 
required. This quantity is usually determined by measurement 
with an empirical curve fit made to the data. Here we shall 
concentrate on boundary layer pressure fluctuations although 
the results will be qualitatively representative of other 
random fluid motion such as engine noise, etc. 

Let us begin by considering the correlation function (16). 
Various authors have given empirical equations. Typical is 
the one given by Y. L. Lin^ 

R p (t; x,y,x*,y*) = V 2 e' a l x *“ x l 

. e -3 | (x*-x)-u c t| e -y I y*-y I C 68 ) 

V 2 , U c , a, 6, y 

are constants chosen to fit the experimental data. Note that 
depends only on the difference of the spatial coordinates. 

To the extent that this is an accurate approximation this greatly 
reduces the amount of experimental information required. 

Let us consider for a moment physical interpretations of 


the constants. From (68) 
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R p (t; x, y,x*,y*) = R p (x;x*-x,y*-y) 


(69) 


For 


x-x* = y-y* = 0, 


R p (t; 0,0) = p 2 


(70) 


Hence, p 2 is is the mean square pressure, see (16), at any 
point on the plate. The mean square at all points is the same 
by virtue of the assumption that only the spatial differences 
appear in (68). 

Next consider |x*-x| ^ 0. The second exponential factor 
of (68) represents a convecting pressure form traveling at 
velocity U c - For fixed x*-x and y*-y, R p will have a maximum 

at T ma* = Cx*-x)/ . 

c 

Finally, the constants a, 0, y may be written 


a = c 1; 

' 6 * 

e = c 2/ * 

Y = c 3/ * 

where 6* is the boundary layer thickness and cj, c 2 , c 3 are 
non-dimensional constants on the order of .01-1. Of course, it 
is expected on physical grounds that the spatial scale would be 
roughly of the size of 6 *. 
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Using the definition of power spectra (18) and the empirical 
expression (68) we have 

00 

$ (oj;x*“X, y*-y) = i / R (T;x*-x,y*-y)e ltJjT dr (72) 

" TT -«> ^ 

= P r 2BU C 1 e - a lx*-x|-Y|y*-y|-i»(x^x) 

7 _(irf 2 +w 2 J c 

This equation is the input needed for (24) or (25) to carry 
our analysis further. Alternatively may be available directly 

from experimental data. 

It should be noted that more complicated correlation and power 
spectra functions are sometimes used. See (36) and Reference [5] . 
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6. CONCLUDING REMARKS 

A review of power spectral methods for determining linear 
response of structures to random pressure fluctuations has been 
given. Various simplifying assumptions are made for the purpose 
of obtaining useful explicit formulae for structural response. 
The transmission of sound through a flexible structure into 
an interior cavity is also briefly treated. 
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